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We obtain the degree distribution for a class of growing network models on flat and curved spaces. 
These models evolve by preferential attachment weighted by a function of the distance between 
nodes. The degree distribution of these models is similar to the one of the fitness model of Bianconi 
and Barabasi, with a fitness distribution dependent on the metric and the density of nodes. We show 
that curvature singularities in these spaces can give rise to asymptotic Bose-Einstein condensation, 
but transient condensation can be observed also in smooth hyperbolic spaces with strong curvature. 
We provide numerical results for spaces of constant curvature (sphere, flat and hyperbolic space) 
and we discuss the conditions for the breakdown of this approach and the critical points of the 
transition to distance-dominated attachment. Finally we discuss the distribution of link lengths. 
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I. INTRODUCTION 

Scale-free networks have attracted a wide interest as 
models for many systems. Their degree distribution can 
be explained by the mechanism of preferential attach- 
ment in growing networks p.]. Preferential attachment 
has been proposed as a realistic effective mechanism of 
network growth for many real- world examples, from In- 
ternet and the World Wide Web to Wikipedia and cita- 
tion networks. However, preferential attachment alone 
is not enough to explain the structure of scale-free net- 
works, since many of these networks build also on other 
variables like fitness 0, 0] and spatial structure J3-[ll| . 

The fitness model of Bianconi and Barabasi [3] is the 
basic example of a network model where the attachment 
depends not only on the number of links but also on 
another variable, namely the quality of the nodes. This 
model predicts multiple power-law scaling in the degree 
distribution and a phase with Bose-Einstein condensation 
on nodes with high fitness Q and it has been applied to 
the WWW iSi] J where Google could be an interesting 
example of an emerging condensate. 

Space is another feature that plays an important role in 
many real-world networks. Examples include the world 
airport network, which is scale-free [l^ and where dis- 
tances are relevant [l^j , but physical and abstract spaces 
appear in many other transportation networks as well 
as in communication and social networks. Another bi- 
ological example is given by brain connectivity in the 
hippocampus, which depends on the distance in phy sical 
space and shows a scale-free degree distribution [16] and 
a short axon length compared to the size of the region. 

It is possible to obtain scale- free networks from models 
based on on hidden variables even without any preferen- 
tial attachment [l^ . Interestingly, there are several 
models with power-law degree distributions induced by 
the space itself (see [l^ for an extensive review) . In some 
models, it is the interplay /competition between spatial 
and network rules that generates to a scale-free distri- 
bution [20l - [2i |. Other models rely only on the spatial 



structure; as an example, recent models based on static 
random networks on hidden hyperbolic spaces [25l - [27| 
have been shown to give rise to scale-free degree distribu- 
tions and could explain the assortativity and clustering 
observed for Internet. 

An alternative, simpler way to obtain spatial scale- 
free networks is to embed the nodes in a metric space 
and modify the preferential attachment rule to includes 
spatial information. The classic preferential attachment 
mechanism is based on a connection probability 11^^^ 
that is linear in the degree of the node ^j^i (x ki. If 
each node i is associated to a position Xi in a metric 
space, an immediate modification of preferential attach- 
ment is given by a multiplicative weight depending on 
the distance dij between nodes, i.e. 



n, 



cx a{dij)ki 



(1) 



The positions Xi can be randomly chosen from a distri- 
bution p{x). If all distance are equal (for example, if the 
space reduces to a single point) then the spatial structure 
disappears and the model reduces to the Barabasi- Albert 
model, otherwise there could be a non-trivial interplay 
between the spatial and network structure. Growing 
network models in this class have been proposed several 
times in the literature and have been studied numerically 
in Euclidean spaces These studies are reviewed 

in JJJ, section IV. D. 2. Furthermore, analytical results 
have been obtained for some models on the sphere [34Ij| 
and in generic spaces under strong approximations [7j, |8| . 
However, despite the simplicity of these models, there are 
no analytical or numerical results for models on generic 
metric spaces and with generic weights a{d). Our aim is 
to provide a general analysis of this class of models. 

In this paper we focus on networks on flat and curved 
spaces, the growth of which follows a preferential attach- 
ment rule weighted by a function of the distance between 
nodes as in equation ([1]). Our main result is the deriva- 
tion of a general expression for the asymptotic degree 
distribution. The form of the degree distribution is es- 
sentially equivalent to the fitness model of Bianconi and 
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Barabasi 0, with a fitness distribution that depends on 
the space and on the connection function cr((i) through 
a selfconsistency equation. We present some examples 
where the distribution can be solved explicitly. These 
examples include homogeneous space and spaces with 
curvature singularities; in the latter, Bose-Einstein con- 
densation can occur on nodes near a singularity. We also 
study numerically the distribution for disks in the flat 
plane, on the sphere and in the hyperbolic plane, showing 
that in the hyperbolic model there is a long but transient 
regime of Bose-Einstein condensation. Finally we give a 
general expression for the distribution of link lengths. 



II. DEGREE DISTRIBUTION OF GROWING 
NETWORKS ON HIDDEN SPACES 



We present a general class of models for networks with 
preferential attachment on a metric space. In these mod- 
els, each node has a position x randomly chosen from a 
space M. with probability p{x). We call Xi and ki the 
position and degree of the ith node. At each time step, 
a node with m links is added to the network and the m 
links are randomly attached to nodes in the network with 
a modified preferential attachment rule. 

We are mostly interested in the case where each node 
is assigned a positions on a Riemannian manifold of di- 
mension D (possibly with boundaries and/or singulari- 
ties) with finite volume, that is, 5* = ds^) where 
is the manifold and ds^ is its metric, that is, the squared 
distance between the points x and x + dx, which is a 
bilinear form ds^ = j gij{x)dxidxj in the infinitesi- 
mal displacement dx on the manifold. Specific models 
in this class were studied in 
general model was defined in 

of [1] is valid only for almost-homogeneous manifolds, as 
discussed in section ITlI Al 

In this work we employ the natural requirement that 
all the information about hidden variables should come 
from the metric structure of the manifold. More specif- 
ically, we require that nodes are assigned random posi- 
tion on the manifold such that the average number of 
nodes per unit volume is constant. Therefore the result- 
ing density is simply p{x) = l/Vol(A^) if expressed in 
terms of (locally) Euclidean coordinates, or more gener- 
ally p{x) — (x) /Yo\{M) where Jix) is the Jacobian 
of the change of coordinates to (locally) Euclidean coor- 
dinates. However, all the results of this section extend 
naturally to the case of a generic node density p{x) arbi- 
trarily chosen. 

Moreover, we require that the preferential attachment 
rule be modified such that the connection probability 
from a new node with coordinate y to a node with co- 
ordinate x depends on their distance d{x, y), that is, the 
attachment rule is given by a real positive function a{d) 
such that the probability to attach a link to the ith node 
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a{d{x,Xi))h 
J2,cridix,Xj))kj 



(2) 



where x corresponds to the position of the new node. 
These models correspond to the heterogenous network 
models of 0, Hi, but while these papers focus on the 
almost-homogeneous case, we allow for all spaces and all 
levels of heterogeneity. 

An important quantity related to cr(d) is its typical 
length scale Xa ~ {da{d)) / {a{d)) , which depends explic- 
itly on the manifold M (and on its boundary, if it exists) . 
Denoting by R the typical length scale of the manifold, 
the connection function is "long-range" ii ^ R and 
"short-range" if Ag- ^ R. The networks associated with 
the two cases can have different properties. 

A network model on a hidden space is therefore defined 
by a manifold A4 and a connection probability (j{d). We 
derive the degree distribution through the rate equation 
approach introduced in [2^, [s^]. The equation for the 
node density (x) , which is the average density of nodes 
with hidden position x and degree k, is 

nk{x,t + 1) = nk{x,t) — m X (3) 
^^^ a{d{x,y)){knk{x,t) ~ (fc - l)nk-i{x,t)) ^^^^ 
M XmC^^^ cr(d(z,y))X;^=™fc'»^fe'(2::0 

plus a term +6k.mp{x) that accounts for the birth of new 
nodes. We assume a linear scaling with time for the quan- 
tity 



d''zaid{z,y)) V k'nk'iz,t) = mCiy)t + oit) (4) 



M 



then can rewrite it as 



nkix,t+l) = nk{x,t)~{knk{x,t)-{k-l)nk-i{x,t)) 



q{x)_ 

t 
(5) 

where g(a:) plays the same role as the (average) fitness of 
the node 0, Q and is defined as 



q{x) 



M 



,D cr{dix.,y)) 
C{y) 



p{v) 



(6) 



The "fitness" q{x) is determined by solving asymptoti- 
cally the above equation 



nk{x,t) = t 



p{x) T{m + q{x)-^)T{k) 
'^T{k + l + q{x)-^)T{m) 

p{x) /fc\-*'^'^"^"^ 



q{x)m \m 



(7) 



and substituting in the definition of C{y) and in ([6|) to 
obtain a single functional equation for q{x): 



q{x) 



d'^y p{y) 



<^{d{x,y)) 



M 



(8) 



3 



This consistency equation determines q(x) as a function 
of the global structure of the manifold S = (A^, ds^), the 
weight function a{d) and the node density p{x). 

From the point of view of the degree distribution, this 
class of models is equivalent to the Bianconi-Barabasi 
fitness model j2j, [S] , but in this case the fitness distribu- 
tion is dynamically determined by p(x) and (j{d) through 
equation ([S]). The resulting degree distribution is 



p{k) 



M 



q{x)m 



(9) 



so the distribution is a sum of power laws similarly to 
the fitness model, and for a regular distribution of x and 
q{x) it typically reduces to a power law with logarithmic 
corrections 

Note that if we allow for an arbitrary node density 
p{x), the only relevant information about the manifold is 
the distance function d{x,y), therefore these results can 
be applied to any metric space. However from now on 
we focus on Riemannian manifolds with uniform node 
density. 

We are particularly interested in warped manifolds. A 
warped manifold is a space M that is locally a product 
RxTi where H is an homogeneous space with metric ds^ , 
and whose metric has the form ds^ = dr^ + f{r)ds'^. 
We are interested in spaces with 0{D) symmetry, which 
means that H is the sphere S^~^ with the natural metric 



dfl . In this case, if there is a radius (say. 



0) such 



that /(r) — > for r — > 0, then /(r) ^ for r ^ to 
avoid singularities in the metric. Flat spaces, hyperbolic 
spaces and hyperspheres can be seen as particular exam- 
ples of these spaces with /(r) — r^, f{r) = sinh^ {(r) / 
and /(r) = sin^(Cr)/C^ respectively. In these spaces the 
0{D) symmetry implies that q and C depend only on r 
and that q{r) is determined through the equation 



q{r) 



where h(r, s) is 



ds 



du 



h{s,u) 



(10) 



hir,s)=f{sy 



dn a{d{r, no; s,n)) (11) 



In the following sections we discuss three analyti- 
cally solvable cases (homogeneous spaces, balls in flat or 
hyperbolic spaces with long-range connections, warped 
throaths) and then the interesting cases of disks on the 
flat plane, the sphere and the hyperbolic plane, which 
are the simplest spaces corresponding to constant zero, 
positive and negative curvature respectively. 



III. SOLVABLE MODELS 
A. Homogeneous spaces 

We consider flrst the cases where Ai is an homoge- 
neous space. This means that there is a symmetry group 



on A4 that leaves the metric (and the measure) invariant 
and moreover that the orbit of each point in M under the 
action of the group is M itself. The first consequence is 
that the volumes d^xp{x) and the connection probability 
a{d{x, y)) are invariant under the action of the symmetry 
group, as are all quantities build from integrals of func- 
tions of the position over the whole manifold M.. Then, 
if the assumption @ is correct and if all the integrals 
involved in the definition of C{x) and q{x) are conver- 
gent, neither C{x) nor q{x) depend on the position x, 
and therefore integrating both sides of the equation ^ 
by / d^xp{x) we obtain q — 1/2. This means that the 
degree distribution of the model on homogeneous spaces 
remains essentially the same as the one of the Barabasi- 
Albert model 



p{k) cx k 



(12) 



For homogeneous spaces, this result does not depend at 
all on the form of (j{d), as long as its integral over the 
volume is convergent. 

The result ([T^ was obtained numerically by 
Barthelemy f3| for the case of the flat disk in D = 2 
with a{d) = e^'^l"^" and Tc not too small. In the same 
paper it was shown that for small the transient regime 
dominates the dynamics for long times, inducing a cutoff 
^raax ~ n^'^ where n is the number of nodes inside a vol- 
ume of order r\. In this section we have shown that the 
result p2)) is actually general for homogeneous spaces. 

Recently, Jordan [ll| proved a slightly more general 
result: the degree distribution (IT^ is the same for all 
spaces where the volume of a ball depends only on its 
radius and not on its position, provided that a{(J) has 
a positive minimum and that the integrals of o(d) and 
<j{d)^ over the whole manifold are finite. This result ap- 
plies to all homogeneous spaces, confirming the above 
argument. Note that Jordan's result can be directly ob- 
tained from equation ([5]) by substituting q{x) = 1/2 and 
interpreting all integrals as Lebesgue integrals, with the 
only requirement d^x (j{d{x,y)) < oo. 

In practice, equation (jl2p is also true for open balls in 
homogeneous spaces if Act <C R, where R is the radius of 
the ball. Consider as an example the hyperbolic disk: if 
the typical radius Xa- of the function a{d) is much smaller 
than the radius of the disk, then the disk appears to each 
node as the (infinite homogeneous) hyperbolic plane. De- 
viations from q(r) = 1/2 decay exponentially with the 
distance from the boundary of the ball. More generally, 
a rough approximation for almost-homogeneous spaces is 
to consider C independent on the position on the mani- 
fold and therefore obtain the approximate equation Q 



q{x) 



2 !m Jj^ dDz p{w)p{z)a{d{w, z)) 



(13) 
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B. Balls with long-range connections 

We consider balls in flat spaces of high dimension D ^ 
1. In this limit, most of the nodes lie within a distance 
R/D from the boundary of the ball. This means that 
for D oo almost all the nodes lie approximately on 
the boundary, that is, on the {D — l)-dimensional sphere 
S^~^, which is an homogeneous space with the induced 
metric. Moreover, if the connection probability grows 
with the distance, most of the nodes will connect to nodes 
on the boundary. Therefore, these nodes have q = 1/2 
and it is possible to calculate q{r) for the few nodes in 
the inside from equation (1101) by considering an effective 
distribution p{r) cx 6{r — R) of nodes. If we choose for 
example a{d) = dP the result is 

/J" M sin(i?)^-2(^2 + ^2 _ 2rRcos{d)) _ R^ + 
^^""^ ~ J^di^ sin(i?)^-24i?2(l -cos(i$l)) " 4i?2 

(14) 

Similar (but more complicated) results can be obtained 
in the same limit 13 3> 1 with other functions (j{d) or with 
balls in other spaces (hyperspheres, hyperbolic spaces). 

C. Warped throats, condensation and the fitness 
model 

A natural question about these models is: what hap- 
pens when equation ([5]) does not admit solutions? The 
absence of solutions is a sign of the breakdown of the 
mean- field approximation in equations (l3]),(|4]). This 
breakdown could appear as condensation of the links on 
a single node as in the fitness model of Bianconi and 
Barabasi 0]. In order to show that this possibility is 
realized at least in some cases, we naturally embed the 
condensate phase of the fitness model in a network on a 
warped geometry. 

The basic mechanism for condensation can be under- 
stood by considering the throat-like geometry illustrated 
in figure [TJ Most of the nodes would lie on the upper 
border of the throat, while only a few nodes would lie 
near the tip. If the connection function cr(d) increases 
strongly with the distance, many nodes would connect 
to the node that is closest to the tip, because it would 
have maximum distance from the border and almost no 
competitors. In turn, the degree of this node would grow 
fast and attract links because of preferential attachment. 
Then, links would condensate on this node. 

For a concrete realization of the mapping between our 
model and the fitness model, we consider networks on 
warped throats. These geometries are well-known in the 
context of string theory and particularly of gauge/gravity 
correspondence [3l|. They can be realized as warped 
spaces x S^^^ with a metric ds"^ = dr^ + f{r)d^'^ 
and a function /(r) that remains very small for a large 
interval of r, as illustrated in Figure [1] In particular, we 
consider r G [0,1], D ^ 1 and we choose /(r) = er*"' 
and a{d) = d. Then p{r) = r'^^^-^^^ikD - k + 2)/2. 




FIG. 1. Example of a two-dimensional warped throat embed- 
ded in tridimensional space. 



In the limit £ <C 1, kD ^ 1 almost all the nodes lie 
at r ~ 1. Moreover, the size of the transverse spheri- 
cal sections is given by ^/e, which is practically negligi- 
ble, therefore the connection probability can be approx- 
imated as a{di^2) — di^2 — \fi — r2\- From the point 
of view of the degree distribution, it can be further ap- 
proximated as cr(c?i_2) — I — ri since r2 would almost 
surely lie around 1. If we define 77 = 1 — r, we see that 
(7(^1^2) — Vi ^-iid therefore the model corresponds pre- 
cisely to the Bianconi-Barabasi model with fitness dis- 
tribution p{r]) oc (1 — 7y)'=(-°-i)/2^ that is, in the deep 
condensate phase [3|. Then the mean- field approxima- 
tion breaks down and there is a condensation of links 
(that is, a finite fraction of the links) on a node near the 
tip of the throat. 

The correspondence between the two models captures 
well the competition between nodes with higher "fitness" , 
that is, the nodes between r = and r' <C 1. This 
means that the condensate is not a transient effect but 
a long-term feature of this geometry. This is related to 
the presence of a metric singularity on the tip of the 
throat: in fact a non-singular tip would look locally like 
E^, therefore it would contain a finite density of nodes 
increasing linearly in time and competing asymptotically 
with the condensate. The competition of nodes with the 
same fitness would make the degree of the most connected 
node grow sublinearly {q < 1) and then the condensate 
would disappear. 

In the above model with fc = 2 it is possible to see 
explicitly the condensation phase transition by varying 
the parameter e . This interpolation is possible because 
in this case e has also an easy geometric interpretation: in 
fact the geometry reduces to a Z)-dimensional hypercone 
with a conical singularity at the tip r = with cone 
angle a and the parameter £ is a function e = /Att'^ 
of this angle. For £ = 1 we have a = 27r, therefore the 
space is a flat ball without singularities and the model 
reduces to the one of the previous section IIII Bl with a 
degree distribution p{k) ~ k~^. On the other extreme, 
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for e — )■ the cone is extremely shallow (a — >■ 0) and we 
recover the model above, which is in the deep condensate 
phase. From geometrical reasoning, the critical point for 
this phase transition in the limit D ^ oo lies around 
Sc ^ 16arcsin2(l/4)/7r2 ~ 0.1. 

It is not clear if there can be condensation without 
a singularity in the metric. For smooth manifolds, the 
manifold looks locally like the hyperspace M.^ in the in- 
side and has an uniform density of nodes everywhere, 
therefore the argument above forbids the appearance of 
a condensate at any position except the boundary. A sim- 
ilar argument should apply to points on a smooth bound- 
ary. For this reason, we conjecture that the equation (|8]) 
for q{x) always admit a solution for a smooth manifold 
of finite diameter with smooth boundaries and a smooth 
function a{d). The argument can be easily extended to 
the case of a generic (non-uniform) density function p{x), 
as long as p{x) is stricty positive and smooth in all points 
of the manifold. 

Note that condensation with cr((i) an increasing (de- 
creasing) function of the distance d occurs typically in 
manifolds with singularities of positive (negative) curva- 
ture. An example of the latter is the above hypercone 
model with e ^ 1 (a strongly hyperbolic space, the op- 
posite limit with respect to a warped throat) and a de- 
creasing function like a{d) = e^'^. 



IV. NUMERICAL RESULTS 

Since it is difficult to solve analytically the equation 
(flUl) . we simulated the above model on some spaces and 
solved numerically the equations for q{x) to compare it 
with the result of the simulations. Simulations were per- 
formed for two-dimensional warped manifolds with con- 
stant curvature: the fiat disk (curvature ( — 0), the el- 
liptic disk (disk on the sphere, C = +1) and the hyper- 
bolic disk (disk in the hyperbolic plane, C = —1). These 
choices of ( are completely general since a rescaling of the 
curvature is equivalent to a rescaling of the radius rmax 
of the disk and of the typical length Xa of the connection 
probability. 

We built the initial network with uq = 2m nodes with 
m links each, randomly connected. To every node, both 
initial and new ones, were assigned random coordinates 
in the manifold with a uniform distribution with respect 
to the volume of the manifold. The rest of the network 
was built as in the model presented above. Data reported 
are for simulation with total number on nodes N = 2 ■ 
10^ and are averaged over 50 simulations with differents 
starting seeds. 

Numerical solution of ([TUl) is based on the numerical 
solution of the system of equations 



ds 



h{r, s)m 



jjiTiax ^^^^ u)k{u, t) 



(15) 
(16) 
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FIG. 2. Plots of q{r) and p{k) (relative to BA) for differents 
a(d) for the flat disk of radius 1.5. 



which is equivalent to a numerical integration of the net- 
work using macronodes. These equations are integrated 
up to a time t when the solution satisfied the condition 



qi (r, i) ~ 1 - m/k{r, t) ^ q2 (r, t) 



(17) 



with a given precision. If the above condition is met, the 
result q{r) is a numerical solution of (jlOp . 

For the simulations, q{r) can be estimated in two ways: 



log 



q{r) 



fci(to)/^._. 



log{tf) - log(to) 
m 



-1 



(18) 
(19) 



where (. . .) j denotes the average over nodes in the set 
/. The difference between q(r) and 92 (^) (or qi{r,t) and 
92 (f", t)) suggests how far the network is from equilibrium. 
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FIG. 3. Comparison between predicted and simulated g(r) 
for flat disk of radius 1 and some choices of a{d). 



A. The flat and elliptic disk 

The flat space is an homogeneous space, so the flat 
disk looks almost homogeneous for a short-range connec- 
tion function a{d) — e"^^'^/^ and q{x) ~ 0.5 except for 
the region near the border, therefore the degree distribu- 
tion of the network is similar to the one generated by the 
Barabasi- Albert (BA) model. Functions like e"'^''/^ and 
especially 9{R—d) have higher A^-, therefore produce net- 
works with a q{x) and consequently a pik) that deviate 
from the BA model, as it can be seen in Figure [2] There 
is a very good agreement between the numerical solution 
of equation ([TU]) and q{r) extracted from the simulations, 
as shown in Figure [3l 

Simulation results for the disk on the sphere are similar 
to the flat case, but the deviations from the BA case are 
slightly less pronounced. We compare the distributions 
of q{r) and p{k) for networks with ^[d) — 0{r„iax — d) on 
disks in flat space, sphere and hyperbolic plane in Figure 

SI 

The case of networks in one- and two-dimensional flat 
space was actually studied in a series of works d, d, [H, 
|32| with the function a{d) = d" for negative values of a. 
This case is interesting because as emphasized in 0, Q ; 
the degree distribution resembles a power-law for a above 
a critical value ac, while below this value it develops a 
stretched exponential tail. This transition to exponential 
behaviour is confirmed by our simulations for the flat disk 
in D = 1 and 2, as shown in Figure [SJ The critical point 
appears to be ac = —1 for D — 1 and ac — —2 for 
D = 2, different from the values found in 4] for the two- 
dimensional torus but in agreement with later results in 
[33] and the theoretical suggestions in 

Our approach offers a simple explanation for this be- 
haviour. For a < —D, the connection function a{d) = d" 
has a divergent integral on a manifold with finite volume 
in dimension D. The rate equation approach leading to 
equation ([S]) breaks down if the integrals involving the 
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FIG. 4. Plots of q{r) and p{k) (relative to BA) for disks with 
rmax ~ 1.5 and (j{d) — 9(rmax — d) in two-dimensional spaces 
of different curvature. 



connection function a{d) diverge. In particular, if the di- 
vergence is localized at d = 0, new nodes connect mostly 
to the nearest nodes and the distribution develops a ex- 
ponential tail. The distribution resembles a stretched 
exponential, as shown numerically for D — 2 in [33| and 
also analytically in a similar model where each new node 
i connects to the existing node with maximum kjdfj f34|. 
Therefore the predicted critical point is ac = —D, in 
agreement with simulations in Figure [5j As an interest- 
ing observation, note that the behaviour of in Figure 
[S] depends apparently only on the ratio a/D and not on 
the dimension D itself, at least for small D. 

Incidentally, the good agreement between the predic- 
tion ac = —D and the simulations shown in Figure [5] 
suggests that all networks in homogeneous spaces sat- 
isfying the condition Jj^d^x a{d{x,y)) < oo have the 
same degree distribution as the BA model and that the 
condition d^x a{d{x, y))"^ < oo assumed in TT] is not 
necessary. 
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B. The hyperbolic disk 

For disks in the hyperbolic plane of curvature -1, the 
results depend strongly on the radius of the disk. The 
deviation from the homogeneous case (BA distribution) 
is more pronounced as Ao- and r^ax (i-e., curvature) in- 
crease. For a disk of radius 1.5, the simulated q{r) and 
p{k) show a pattern slightly more pronounced than the 
flat case, while results for a disk of radius 4 are shown in 
Figure [B] 

The hyperbolic space is quite interesting because it 
resembles the scenario opposite to the warped throath 
discussed above, given the exponential expansion of the 
volumes as r increase, therefore it is a good place to ob- 
serve condensation in models with decreasing cr(d). Ac- 
tually, the absence of singularies forbids the existence 
of an asymptotic condensate, but transient condensates 
with q « 1 could last for exponentially long times due to 
the exponentially small number of competing nodes near 
the center of the disk. This effect can be clearly seen in 
Figure El with a{d) = 9{rmax - d). 



The fact that q{r) goes above 1 (that is, superlin- 
ear scaling) depends only on the long out-of-equilibrium 
transient regime, as it can be seen by comparing all esti- 
mators of q(r) in Figure[7l For large times, g(0) is slightly 
below 1. 



V. DISTRIBUTION OF LINK LENGTHS 

The asymptotic distribution of link lengths follows the 
equation 

vil) = f d^x f y))p{x)p{y)5{d{x, v) - 

(20) 

which depends on the fitness q{x) defined by equation 
For (almost) homogeneous spaces, the above equation 
reduces to the simple distribution 



P{1) 



<l)Pd{l) 



j^dVa{V)p,{V) 



(21) 
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FIG. 7. Plots of q{r) calculated numerically (both at finite 
and asymptotic times) and simulated for the hyperbolic disk 
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FIG. 8. Link length distributions p{l) averaged over 50 simu- 
lations for the flat disk of radius 1 and connection probabilities 
a{d) = 6(1 — d) and e~^^'^. The continuous lines represent the 
predictions p{l) = 216(1 — I) and 256/e~^^' for homogeneous 
spaces. 



where pd{l) is the density of nodes at distance I from 
a random node. This result is consistent with previous 
studies with a{d) = d" in Q- 

Two examples of link length distributions in flat space 
are shown in Figure [S] The model with short-range in- 



teractions a{d) 



is almost homogeneous and fits 



well the prediction p{l) = 256Ze~^^' from equation ([21]) . 
while the one with long-range interactions a{d) = 9(1 — d) 
deviates from the distribution p{l) — 216(1 — I) expected 
for an homogeneous space. 



VI. CONCLUSIONS 

The results of this paper show that the power-law be- 
haviour of the degree distribution of growing networks 
with preferential attachment is quite robust with respect 
to the spatial structure of the network. The effect of 
the interplay of preferential attachment and the spatial 
nature of the network results in an heterogeneity be- 
tween nodes in different positions. In spaces with high 
curvature or singularities, it is possible to observe Bose- 
Einstein condensation driven by the intrinsic heterogene- 
ity of the space or by the finite size of the network, as it 
happens in hyperbolic spaces at strong curvature. 

These results have been obtained through a rate equa- 
tion approach, which breaks down if the integrals in- 
volved in the selfconsistency equation for q(x) diverge. 
This is the case for connection functions a(d) with di- 
vergent integral. In the simple case (T(d) = d", there 
is a critical point Uc = —D for the transition from the 
scale-free behaviour to an exponential behaviour due to 
distance-dominated attachment, that is, preferential con- 
nection to the nearest neighbours. 

In this paper we only considered a spatial embedding 
of the basic Barabasi- Albert model, which gives rise to 
a power-law degree distribution with exponent 7 = 3. 
However, most scale-free networks found in natural and 
technological systems have a degree distribution with an 
exponent 7 < 3. There are several variations on the 
(see 



for a review) that are 



Barabasi- Albert model 
able to account for a different exponent and that can be 
included naturally in the spatial models studied in this 
paper without changing the main results. For example, 
adding m' extra links per unit time would leave equa- 
tion ([9|) invariant but change the l.h.s. of equation (|8]) 
to q(x)/{l + 2m,' /m), which corresponds to an exponent 
7 = 2-1- 1/(1 -l- 2m' /m) in homogeneous spaces, varying 
between 2 < 7 < 3. 

Finally, spatial networks have other interesting proper- 
ties, like clustering and assortativity, that have not been 
studied in this paper and would deserve further work. 
In the models considered here, the volume Vol(A^) of 
the manifold is constant and the density of nodes grows 
linearly in time, going to infinity in the thermodynamic 
limit. This divergence in the density makes the cluster- 
ing decrease in the thermodynamic limit. An interesting 
variation on these models would be a growing network 
with constant node density on an expanding space. 
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